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Abstract 

We estimate the stopping length of an energetic gluon in a thermal plasma of strongly 
coupled M = 4 super- Yang-Mills theory by representing the gluon as a doubled string rising 
up out of the horizon. 
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1 Introduction and summary 

Following [TJ El [3], Hj, considerable effort has been devoted to understanding what string 
theory might have to say about the physics of hard probes in heavy-ion collisions. The 
key question is how quickly a hard parton loses energy as it passes through the quark- 
gluon plasma (QGP). The approach of [1] (see also [5]) is tied to the BDMPS jet-quenching 
formalism j6j U\ [8] (see also P, [10], and [H] for a recent review), with a definition of q in 
terms of a partially lightlike Wilson loop. The approaches of [H [31 H] are limited to heavy 
quarks and focus on drag and stochastic forces. Here we would like to propose an extension 
of the approach of (21 0] to accommodate gluons in M = 4 super- Yang- Mills theory (SYM). 
An energetic, off-shell gluon in the thermal medium should be represented as a doubled 
string coming up out of the horizon of the AdSs-Schwarzschild geometry. Such a string must 
eventually fall back into the horizon. So the question is how far the string gets before it 
does so. We will estimate the maximum penetration length Ax, as a function of the initial 
energy E, which is assumed to be much greater than the temperature. If we define 

A IE 
x = nix E— — — , (i) 
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where T is the temperature and gy M iV is the 't Hooft coupling of SYM, then the relation 
we find, for sufficiently large E, is 



Ax « 0.95£ 1/3 . (2) 



This relation is obtained from averaging the leading behavior of the analytic estimates ( 1361) 
and (139]) . The scaling Ax oc E 1 ^ 3 is different from the BDMPS scaling Ax oc E 1 ^ 2 , but not 
very different. In order to compare with BDMPS, we make a rough operational definition of 
q in terms of the stopping length Ax of a gluon of energy E: 

= AE 

grough = 3a s (Ax) 2 ' (3) 

Comparing with BDMPS energy loss is hazardous because the underlying physical picture is 
significantly different. Nevertheless, we plug our estimates of Ax into ([3]) to obtain estimates 
for the corresponding value of q. To extract numerical values, we consider QCD at a temper- 
ature of 280 MeV, which is representative of central gold-gold collisions at ^/snn = 200 GeV. 
We exhibit our estimates in figures [6] and These figures differ only in how we compare 
M = 4 SYM to QCD: in the nomenclature of [12], figure [6] uses the "obvious" scheme and 
yields q ~ 92GeV 2 /fm in the range E = 5 — 25 GeV for the gluon; and figure [7J uses the 
"alternative" scheme and yields q ~ 21GeV 2 /fm in the same energy range. For reasons 
explained in [12], we prefer the alternative scheme, where comparisons are made at fixed 
energy density and the coupling is chosen to make the quark-anti-quark potential in M = 4 
SYM agree as well as it can with lattice results for QCD at separations on the order of 
0.25 fm. In any case, it seems clear that energy loss and thermalization as estimated from 
our falling string picture is more rapid than in the BDMPS formalism with q taken either 
from perturbative estimates or from [I]. On the other hand, according to [13], comparison 
of parton quenching model calculations [T4"l 1X5] to PHENIX data leads to the following 3a 
range for the averaged value (q): 

GeV 2 . /A . . GeV 2 m 

7— < (q <28— , (4) 

tm fm 

with lowest x 2 at (q) w 13 GeV 2 /fm. It is pleasant that our estimate of q using the "alterna- 
tive" scheme falls well inside the experimentally favored range (jlj). However, we emphasize 
that there are significant caveats to this comparison, to be discussed further in section [7J 
The idea that an off-shell gluon in a thermal medium should be represented as we have 
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suggested has several antecedents, including [16j EZ] O I n PS] it was argued that in computing 
Wilson loops at finite temperature, the configuration of two anti-parallel strings rising from 
the horizon to the boundary, and representing a widely separated quark-anti-quark pair, 
receives a color factor of N 2 because of the two string ends on the boundary^ This makes 
sense because the N D3-branes are in some sense "behind" the horizon, and a string ending 
on one of N D3-branes indeed acquires a fundamental or anti-fundamental charge. We had 
been in the habit of thinking of color charges living on the boundary of ArfS'5-Schwarzschild, 
but it seems more faithful to the D-brane origin of a near-extremal black 3-brane for the color 
factor to come from ends at the horizon. Seen in this light, the trailing string of [2j H] derives 
its fundamental color charge from the fact that it actually passes through the horizon. (See 
[20] for a particularly clear exposition of the geometry of the trailing string.) What could 
be more natural, then, than to turn the heavy quark into an energetic gluon by letting the 
string double over on itself and pass back down into the horizon, rather than rising all the 
way up to the conformal boundary? 

The work of [17] employs a zero-temperature limit of approximately this construction 
to consider collisions of gluons. In pure AdS$, however, one can insist upon the view that 
color degrees of freedom "live" on the boundary: upon conformal compactification to global 
AdS§, the endpoints of the strings are seen to rise back up to the boundary at a point which 
is infinitely far from the collision region. In [21], some results of [17] were extended to finite 
temperature. Whereas in [T7| |2T] the focus was on scattering amplitudes of several hard 
gluons, here we are interested in the propagation of a single hard gluon through the thermal 
medium. 

The organization of the rest of this paper is as follows. In section [2] we explain how to 
calculate the energy of a gluon represented as a doubled string rising vertically up from the 
horizon. Section [3] shows how to carry out an analogous computation when the shape of 
the string is part of the trailing string. Section H] detours into the computation of lightlike 
geodesies in the Ad^-Schwarzschild geometry. Section [5] presents estimates of Ax as a 
function of E using the lightlike geodesies discussed in section|H Section [6] presents estimates 
of Ax as a function of E using lightlike geodesies on the worldsheet of the trailing string. 

1 Strings falling into anti-de Sitter space have been considered in other contexts related to heavy- ion 
physics in HUGS] . 

2 When the strings are separated, they only interact by virtual exchange of closed string states, and such 
exchanges are suppressed by powers of N . The main point of the counting powers of N in [16] was that the 
N 2 from the horizon cancels against a 1 /N 2 suppression of this type in order to produce a final amplitude 
which is parametrically comparable to the one coming from a string that joins the two quarks without passing 
into the horizon. 
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Section [7] describes comparisons with the BDMPS energy-loss formalism, expanding on the 
brief discussion above. Section [S] describes a lightlike limit of the falling string which is 
analytically tractable. We end in section [9] with some discussion of possible extensions of 
the falling string picture. 



2 Estimating the energy of a doubled string 

If we accept that an off-shell gluon traveling through a thermal medium should be represented 
as a string with both its endpoints passing through the horizon of AdS^-Schwarzschild, the 
next question is what the shape of the string should be. A natural first guess is that it should 
be straight up and down. The problem is that a string that is straight up and down at t — 
will not hold its shape as it moves in the positive x 1 direction. This is illustrated in figure [TJ 
To demonstrate that the string can't stay vertical, recall first the AdS's-Schwarzschild metric: 

L 2 ( dz 2 \ z A 

ds 2 = Gn V dx^dx v = — I —hdt 2 + dx 2 H — — ) where h = 1 T , (5) 

z z \ h 



and the Nambu-Goto action: 

Sng = _ ^~7 J d2(T V-9 where g af3 = Gf, u d a X^d p X u . (6) 

Our convention is to use indices a, (3 for the string worldsheet coordinates a a , and capital 
X M = X^{a a ) for the embedding coordinates of the classical solution under consideration. 
The worldsheet current of spacetime stress-energy is 

^ = —^8^, (7) 

and the equations of motion following from ([6]) are V a Pf? = 0. Here V a is the covariant 
derivative with respect to the worldsheet metric g a p. Five-dimensional indices like \i are 
treated as scalars with respect to V a . Suppose we start the string in a straight up-and-down 
configuration at time t = 0, as illustrated in figured], with an initial velocity profile v = v(z) 
in the x 1 direction. Using coordinates o a = (t, z), this means that, at t — 0, 
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Figure 1: If a string starts at t = in a straight up-and-down configuration, it doesn't hold 
its shape as time evolves forward. At a later time, indicated as t = 1 in the figure, one must 
solve difficult classical equations of motion to find the shape of the string. But because of the 
infinite redshift characteristic of black hole horizons, the point where the string comes out 
of the horizon cannot move at all. At times t > 1, the string continues to fall down toward 
the horizon. Although it takes an infinite time to fall all the way in, it only propagates a 
finite distance Ax forward. 
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Already from ([8]) we can see that there is a depth- dependent limit on the velocity: v < \fh 
in order for the metric to have Lorentzian signature. Variants of this speed limit have been 
discussed in a number of papers, including [21 [221 [231 [2H [23 [2S]- It implies that the string 
indeed cannot stay upright: it would do so only if v is constant in z, and the only way that 
can be true is if v — 0. To put it another way, the point where the string comes out of the 
horizon cannot move in the x l direction. All the string can do is to fall over into the horizon. 

Assuming the speed limit v < \fh is satisfied everywhere along the string, the total 
momentum of the string can be computed as 



2ixa' 



dz 

^2 



h 



I 



(9) 



The first four components p m of can be identified with the four-momentum of the boundary 
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gauge theory. This is because they are momenta defined in reference to Killing vectors d/dx m 
for m = 0,1,2,3. (We use mostly minus signature, so the energy is — p > 0.) The fifth 
component p z does not have such a simple interpretation. The integral in (Q should be taken 
over the intersection of the string worldsheet with the hypersurface t = 0. More specifically, 
it should be taken over that part of the string worldsheet that is outside the horizon. Because 
this part is doubled over and rises to a minimum depth zuy, we find 



There are two subtleties that affect ([TO]) : 

• The kink at z = Z\jy can support a finite lightlike momentum Sp^, which would have 
to be added to p M . Light-like means G^ u 5p^5p u = 0, which is to say h(5p ) 2 = (Spi) 2 
if 5p^ = for n > 1. 

• One could choose to run the integral over the part of the worldsheet behind the horizon. 
Not doing so is a physical choice, motivated by the fact that nothing behind the horizon 
can classically influence what's outside. We regard whatever the string does behind 
the horizon as part of the dynamics of the thermal medium. 

The result ffTOj) is analogous to the expressions E = m / -f and p = mv-f for a massive particle. 
It shows that there are qualitatively different ways in which to make the string representing 
the gluon highly energetic: one may either take zyy — > 0, or make the local "Lorentz" factor 
1/Vh — v 2 big over some portion of the string worldsheet. 

It has recently been emphasized in [27] that a quasi-particle description of the QGP may 
not be valid. A quasi-particle picture is even less likely to capture the physics of strongly 
coupled M = 4 super- Yang-Mills, where the weakly coupled degrees of freedom are manifest 
only in the dual gravitational description. On the other hand, as long as a gluon has energy 
and momentum much greater than the temperature, it scarcely notices the thermal bath, 
and there should be an approximately unique way to describe it. The small z region of 
AdS^-Schwarzschild is where the presence of the horizon doesn't matter, and it is associated 
with UV physics because an object there translates into a tightly localized or highly energetic 
configuration in gauge theory. So we recover the intuition that a hard gluon can be defined 









(10) 
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with little ambiguity by assuming that most of its momentum comes from the small z region: 



z integration: that is, we ignore the bath altogether. Depending on the context, a better 
approximation may be needed: for instance, one might want to replace 1 — v 2 by h\jv — 1,2 
so that the limit on v previously discussed is correctly implemented for finite but small zjjy. 

The result f JTTT) also focuses attention on the fact that the string describes an off-shell 
object: 



Recall that we use mostly plus signature: thus with the explicit sign included in (112j) . Q 2 > 
means timelike momentum. We do not know how to represent a gluon with Q 2 < 0. The 
limit v — > 1 from below with p held fixed corresponds to taking the gluon on-shell. But 
this limit is not available for T > 0, because eventually it would force z\jy to become 
greater than zjj, and ignoring the bath would then be wrong. This discussion highlights 
a crucial difference between our approach and the more conventional BDMPS treatment, 
where the first step is to make an eikonal approximation where the gluon travels strictly at 
the speed of light. While such an approximation is reasonable in perturbation theory and 
makes sense for sufficiently energetic probes of a finite-sized medium in an asymptotically 
free theory like QCD, we harbor some doubts about the consistency of expanding around 
light-like trajectories of charged particles in theories such as M = 4 super- Yang-Mills where 
the coupling is finite even in the ultraviolet. In any case, a particle which propagates only 
a finite distance through the medium is not on-shell because it is not an asymptotic state. 
Thus we are more reassured than concerned over being forced to take v < 1. 

3 Energetic gluons and the trailing string 

With an approximate expression fTTTj) in hand for the momentum of a hard gluon repre- 
sented by a doubled string, the next thing we should ask is how far in the x l direction the 
doubled string travels before it falls into the black hole. Falling into the black hole corre- 
sponds to thermalization in the dual gauge theory. So if the string travels a distance Ax 
before falling in completely, we expect that Ax should be identified, at least approximately, 




The expression comes from setting h = 1 and Zh = oo in fTTDT) before carrying out the 




(12) 
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as the stopping distance of the hard gluon in the medium. Both the energy, — p Q , and the 
penetration length, Ax, are functions of v and zuv> so the question of maximum penetration 
length can be phrased as maximizing Ax with — p held fixed. This extremization problem 
is challenging because evaluating Ax involves solving the non-linear (but classical) equations 
of motion for the string, starting from an initial state which is only approximately specified. 
Instead of tackling this problem head-on, let's go back to the observation that the straight 
up-and-down string configuration considered in section [2] doesn't hold its shape for t > 0. 
Is there is some other string configuration which does? The trailing string of [2j H] suggests 
itself immediately: it is a steady state solution of a string whose endpoint on the conformal 
boundary is required to move with a definite velocity v. That shape is specified in the gauge 
a a = (t, z) by the embedding 

X 1 = v [t + t(z)} where £(z) = flog + i log ±±l) . (13) 

Here we have introduced a rescaled depth coordinate 

y = j-. (14) 

We will persist in using worldsheet coordinates a a = (t,z), but we express results, such as 
the right hand side of ffT31 . in terms of y when convenient. 

Studies in [2] (see also [28]) show that the trailing string is stable against small pertur- 
bations. As an alternative initial condition to the straight up-and-down string discussed 
previously, let's therefore consider a string which starts at t = in the shape (fT3|) . except 
that it rises to a finite minimum depth z\jy before doubling back over itself and going back 
down into the horizon. See figure [2J What we will find is that although this string doesn't 
hold its shape exactly, it "almost" does, in a sense we will describe in section |6l when its 
momentum is large. 

The worldsheet current of spacetime energy-momentum on the trailing is computed as 
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Figure 2: A string starts at t = in the shape of a falling string that extends up to a 
finite minimum depth, y = ?/ uv . (Recall that y = z/z H .) If the string extended up to the 
boundary, as shown with a dashed curve, its endpoint would move at a speed v. The lightly 
shaded region, below the trajectory labeled POND, retains the shape of the trailing string. 
At times t > (for instance at the time denoted t — 1 in the figure), the string probably 
projects somewhat beyond the POND trajectory into the narrow region between it and the 
null spacetime geodesic. The POND trajectory and the spacetime geodesic are mutually 
tangent at the point t = 0, y = y\j\. 



follows: 



dX* 
da a 



(\ 

v 







-vy 2 /h 


1 
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-h + v 1 



—vy/h 



zW \-v 2 y 2 /h (h + v 2 - hv 2 )/h 



2 A„,2\ IU2 



1 



-h-v 2 + hv 2 v 

—hv 2 y 2 hvy 2 —h + v 2 



(15) 



The total momentum of the doubled string may be computed at t = in a fashion analogous 

to m-m-. 



dz y^gP. 



Jl 1 f 1 dy 

7ra'z H y/\ 



1 f d]J_ f_ h _ 2 h 2 v2y2\ 



(16) 



It is important in (TIB"]) that we stipulated that t = 0, because for t > the string will not 
hold its shape precisely: more on this later. 

The four-dimensional components p m of the momentum have a logarithmic divergence at 
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y = 1. (The last component, p$, diverges as a power, 1/(1 —y).) In the original context of a 
heavy quark propagating through the thermal medium [2J, this infrared divergence owes to 
the fact that the trailing string is the shape attained in a late-time limit, and it accounts for 
the large amount of energy that has already been transferred (or mostly transferred) from 
the quark to the thermal bathjfl In the current context, the divergence should be regulated 
somehow, because we have in mind creating a gluon with large but finite energy at t — and 
then asking how far it propagates. The simplest regulator is simply to cut off the integral in 
([IB]) at some yi& slightly less than 1. 

Instead of measuring the energy and momentum of the string at t = 0, one can obtain 
a divergence-free definition of p m by calculating the amount of energy and momentum in 
the string that makes it past a fixed value of x . More precisely, we should compute the 
flux of worldsheet energy-momentum P® through the intersection of the string worldsheet 
with the hypersurface x 1 = constant. To prepare for this computation, consider a general 
conserved worldsheet current Q a . Its flux through a curve S on the worldsheet, specified by 
er a = er a (?7), is 

f da 13 
Q = j^dr)V=ge a/3 Q a —, (17) 

where e a p is the antisymmetric tensor normalized so that €\2 = 1. Now let £ be the curve 
at x 1 = on the trailing string, as shown in figure EJ Using r\ = z as a parameter for this 
curve, and replacing Q a by P™ in (I17jl . one obtains 



p^ Ax = dz^g 



Pt-l-) P° 



zuv 

L 2 f 1 dy 

yw 

dt 



na'zH y/l - v 2 J yvv y 



1 v £ . 



where in the second expression the derivative (§|) i is taken at constant x 1 . Unlike in ( fl6l) . 
the integrals defining the four-dimensional components p^ cd x are convergent, so they don't 
require an IR cutoff. 

We are reassured to observe that for ?/uv -C 1, we recover the form (fTTj) from both 016|) 
and ffTS]) . We will ignore the possibility of an additional lightlike contribution to p^ from the 
tip of the string, where it doubles over. 



3 Thc infrared divergence provides an extreme example of how various shapes of the string describe various 
states of the gluon and the bath. The infrared tail of the trailing string encodes the interaction of the gluon 
with the medium in such a way as to form a sonic boom plus a diffusion wake. 
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4 Lightlike geodesies in AciS'5-Schwarzschild 



The tip of an open string, or a doubled string, must move at the speed of light. But it 
usually does not follow a lightlike geodesic, because it is being pulled in some direction by 
the rest of the string. In the case under consideration, the pull is downward (in the positive 
z direction) and backward (in the negative x 1 direction). So to find an upper bound on how 
far the string gets in the positive x l direction before falling through the horizon, we could 
consider the trajectory of a lightlike particle that starts at the tip of the string at t — and 
falls into the horizon without experiencing the pull of the string. To this end, let's work out 
free particle trajectories in the AtiS's-Schwarzschild geometry. Parameterizing the particle's 
worldline A M = X^{rj) with an arbitrary variable rj, the action may be expressed as 



S = \ I 



1 dX»dX v 2 ' 

-G^v—, -, me 

e arj drj 



(19) 



Here m is the mass (eventually to be taken to 0) and e is a Lagrange multiplier. For m ^ 0, 
one may use the constraint equation for e to eliminate e from the action. After doing so, the 
action ( fl9l) reduces to the standard one, 

S = - I dsm. (20) 



An advantage of (fl9|) is that its m — > limit correctly describes the dynamics of massless 
particles. 

Let's work in a gauge where rj = z and consider trajectories of the form 

X° = X°(z) X l =X\z) X 2 = X 3 = 0. (21) 

Then 

S = [dzC where C = ^— ( -h(X ') 2 + (A 1 ') 2 + t V ~™ 2 e , (22) 
J 2ez 2 \ hj 2 

where primes indicate d/dz. One may immediately form conserved momenta 

» s ^ = -5**° » s £s = 7* xV - (23) 

The equation of motion for e is an algebraic constraint: 

L 2 1 



Z 2 y^2 _ hp 2 _ hm 2 L 2/ z 2 



(24) 
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-1.5 -1. -0.5 0. 

x l /z H 

Figure 3: Typical orbits for a massless particle in AdS^-Schwarzschild, all leading into the 
horizon (the dashed black line) at x l = 0. An open orbit is shown in red; the critical orbit 
is shown in green; and a closed orbit is shown in blue. 



In order to make energy positive, po should be negative, so we should choose the plus sign 
in ( |24l) for trajectories (or segments of trajectories) where z increases as t increases, and the 
minus sign for trajectories where z decreases as t increases. Hereafter we will always choose 
the plus sign, corresponding to particles falling down toward the horizon. 

The particle trajectories X M = X^(z) can be determined using (1231) - (I2ll) once one spec- 
ifies po and p\. In the massless limit, the shape of the orbits X 1 = X 1 (z) is determined 
from 

dX 1 u ez 2 Pi/pq 
dz L 2 v 71 " h v\lvl 

To perform the integral in (1251) one needs elliptic functions, and the explicit result is not 
very enlightening. But for pi = —po, the result is very simple: 

^L = K-^- = K--, (26) 
z H z y 

where K is a constant of integration. One may also straightforwardly show that 

X^X 1 -^), (27) 

where £(z) is as defined in (Tl3|) . We will refer to (126]) as the critical orbit, because for pi < —p 
the orbits intersect both the conformal boundary and the horizon, while for pi > —p the 
orbits begin and end at the horizon: see figure In all cases, getting to or from the horizon 
takes an infinite amount of coordinate time X°. 
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5 Estimating the penetration depth using spacetime 
geodesies 

To find an approximate upper bound on the penetration depth, we should start a massless 
particle at the tip of the string with a physically motivated choice of p^. But how should we 
make this choice? To answer this question, let's have a look at the geometry of the trailing 
string in the hyperplane x 2 = x 3 = 0. Because the intrinsic geometry is Lorentzian, there is 
a basis of lightlike vectors, £ a and k a , on the worldsheet: using coordinates a a = (t, z) as 
before. 



h(h-v 2 ) 
1 



k ° = I I • (28) 



All components of i a and k a are non-singular functions of z, and all are positive except for k z 



which is positive only when h < v 2 , meaning below the worldsheet horizon at z — ZhvT—v 2 . 
One may correspondingly form spacetime vectors = i a dX^ / 'da a and = k a dX^ / 'da a . 
To complete a basis for the x 2 = x 3 = hyperplane one may add the vector n^ 1 normal to 
the worldsheet: it is defined up to an overall factor by the equations n^dX^ /da a = 0, or 
equivalently = = 0, and one easily finds 

n M = (-v 1 . (29) 

The vector points in the direction of a lightlike signal traveling down the string. A sensible 
initial condition for the massless particle whose trajectory is supposed to approximately 
bound the motion of the string worldsheet is p^ oc i^. Combining this initial condition 
with fl25|) . we can immediately calculate how far the massless particle gets in the positive x 1 
direction before falling into the horizon: 

Ax spacEtime = -z H f dy- _ where 

Jwv V 1 - (! -y 4 )pi/Po , 3Q . 

Pi _ V 7 ! - v 2 - y%y 



Po v^lv - (1 - y* v )Vl - v 2 

As remarked previously, the integral can be done in terms of elliptic functions, but the explicit 
form is unenlightening. The subscript "spacetime" in (1301) reminds us that the calculation 
hinges on lightlike geodesies in the ArfS'5-Schwarzschild spacetime. 

The extremization problem that we set out to solve was to maximize the distance Ax 
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traveled by a classical string with fixed energy E = —po, starting from an initial configuration 
with both ends passing through the horizon. What we can now do instead is to maximize 
Ax S p ace ti m e subject to fixed energy. We will consider three ways of defining the energy: 



which comes from ffTTj) ; 



E 



trailing 



ira'z H vT 

which comes from ffl~6l); and 



T 2 1 1 
Evv = — 7^^— , (31) 



-^fixcd x 1 = — ; 75 = ( 1 J ' ( 33 ) 

na'z H VI - v * \yuv J 

which comes from (fTHI) . Evidently, the dimensionful parameters L, a', and zjj enter into 
these expressions only as multiplicative prefactors. So it is convenient to scale them out as 
in (CQ) by defining 

~i x 1 m a 7ia'z H ^ IE 

(34 ) 

where we have used the standard relations 

1 L 2 

Ztf = — = 3ym^ (35) 

Til a 

for A/" = 4 super- Yang- Mills theory. The dimensionless quantities A£ space time, ©uv, and 
-S'fixed x 1 are functions only of v and y\jy, and -^trailing is a function only of v, y\jy, and ?/t R . 
Let's regard i/tr as a fixed cutoff. Then the extremization of A£ spacetime with either Ejjy, 
-^trailing, or -E fixed x i held equal to some fixed value E is a well-defined problem, and we 
shall denote the result Axa{E). The index A labels the assumptions that went into the 
calculation. For example, if we held E\jy fixed in an extremization of Ai spacct i me , then we 
would say A = {spacetime, UV}. If instead we held -©trailing fixed, say with y m = 0.9, then we 
would say A = {spacetime, trailing, ym = 0.9}. In figure H] we show a number of evaluations 
of Axa(E) for several different choices of assumptions. 

As can be seen from figure HI the functional dependence Axa{E) is similar for the various 
sets A of assumptions. To better understand this functional dependence, we computed a 
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■ {spacetime, UV} 

▲ {spacetime, trailing, yiR=0.9) 

• {spacetime, fixed x 1 } 

• {POND, trailing, y m =0.9, cutoff} 

• {POND, fixed x 1 } 



Figure 4: Evaluations of the maximum penetration depth Axa{E) with a variety of assump- 
tions. The solid blue curve shows the analytic approximation ( 139]) to the blue circles, and 
the solid black curve shows the analytic approximation (1361) to the black diamonds. 

series expansion at large E for A = {spacetime, fixed a; 1 }: 

Ax spacctimCifixcd x i = 1.0185E 1 / 3 - 0.8180 + 0.052E- 1 / 3 + 0.017E" 2 / 3 + O^" 1 ) . (36) 

As shown in figure HI this expression provides a fairly good approximation to the numerical 
evaluations. We defer the explanation of how we derived (I3"61) until the end of section [HI 
where the same method will be described in a slightly simpler setting. 

To illustrate more explicitly the nature of the extremization problem, in figure [5] we plot 
Axspacetimc as a function of 7 = 1/yl — v 2 at fixed E &xed x i. We also plot Axpond, a quantity 
to be explained in section [61 It is worth noting that the maximum in Arr is fairly broad in 
7, and it gets broader as E increases. 

6 Estimating the penetration depth using worldsheet 
geodesies 

As we have described, a lightlike geodesic emanating from the tip of the doubled string 
provides an approximate upper bound on how far its trajectory reaches forward in the x l 
direction — approximate because the selection of initial direction is a physical choice based 
on our expectations of what a typical initial state of the string should be. There remains 
the possibility that the string will always fall into the horizon much more quickly than the 
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Ax 




Figure 5: The penetration length Ax as a function of 7 for fixed E. The blue curves represent 
the penetration length for a spacetime geodesic, while the black curve represent the same 
quantity for the point of no disturbance (POND), which is explained in section [6j For all 
curves, the energy was computed at fixed x l using (1331) . The red points mark the maxima 
of each curve and correspond to data points in figure HI 
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lightlike geodesies do. Ideally we would like to have a lower bound to show that this does 
not happen. But recall the nature of the extremization problem: the maximum penetration 
depth Ax(E) is the furthest that a string can get with specified initial energy. It's easy to see 
how it could get much less far: by taking both v and y\jy very small, with total energy held 
fixed, one finds configurations that fall almost straight down into the horizon despite having 
a large energy. What we really want, then, is a demonstration that there is some string 
configuration which gets almost as far as the lightlike geodesies we have already studied. 
This turns out to be straightforward, as we will now discuss. 

The configuration we want to study is again the doubled trailing string shown in figure [2j 
We first mentioned it in connection with the hope that it would almost hold its shape as it 
falls into the horizon. Here's an argument that it does. First imagine the original trailing 
string of [21 0] , which continues all the way up to the conformal boundary. By construction, 
it holds its shape for all time. Now, at time t = 0, let's cut the string at some finite z\jy, 
and let's keep track only of the lower part of it (i.e. the part with z > z^jy). For some time t 
slightly greater than 0, most of this lower part of the string must be in the same shape that 
it would have been had we not cut the string, simply because it hasn't "figured out" that 
the cut occurred. More precisely, if a lightlike signal hasn't had time to propagate down the 
string below a certain depth z*, then the part of the string with z > z* must be in the same 
configuration that it would have been had we not cut the string. The motion of the part of 
the trailing string which the lightlike signal along the worldsheet is able to reach could be 
complicated, and we will not try to figure out what it is, except to restate that it should be 
bounded by the lightlike spacetime geodesic studied in section A lightlike geodesic on the 
string worldsheet, not in bulk spacetime, is relevant for finding out which parts of the trailing 
string remain undeformed at a given time t > because particles like gravitons traveling on 
lightlike spacetime geodesies can affect the string only through string interactions, which are 
suppressed by 1/N. 

We already calculated the tangent vectors to lightlike geodesies in f )28|) . The relevant one 
is £ a because it goes down the string worldsheet (or goes down "faster" in the region where 
h < v 2 ). The differential equation dt/dz = t l is readily integrated, and it gives 

— = t log + - log : y — log r-Z— , 37 

z H 4 l-y 4 1-iy 2 1- 

where 7 = 1/^/1 — v 2 as usual. Plugging (1371) into (113j) . one finds the orbit x 1 = x 1 (z) in 
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spacetime of the point below which the trailing string must be undisturbed: 



x 



POND 



z H 2i 1 - iy,fy 



POND stands for point of no disturbance. A^pond is the difference between x PO nd when 
y = 1 and its initial value when y = yuy. 

Whatever the string may do above the point of no disturbance, Axpqnd sets a lower 
bound on how far it gets in the x l direction. As is evident from figure [51 Axpond is almost 
as big as A£ spacetimc for v close to 1. And because Axpond is a function of v and yuy, it can 
be passed through the same maximization procedure as described at the end of section [5] 
to give estimates Axa(E) of the maximum penetration depth, where A includes "POND." 
A caveat is that A£pond involves an evaluation of the point of no disturbance in the limit 
y — > 1, where it crosses the horizon; but the estimate trailing of the energy of the string 
requires an infrared cutoff at some finite value yip < 1. Evidently, there is an inconsistency 
in whether or not we include in the calculations the part of the trailing string closer to the 
horizon than yip. This caveat can be avoided by using the fixed x l strategy for evaluating the 
initial energy. But to check that it is unlikely to influence our qualitative conclusions when 
we use an infrared regulator to calculate the energy, we considered an alternative definition 
of A£pond which is the difference of £p OND evaluated at ?/ir and ?/uv rather than at 1 and 
yuv- This is indicated by including "cutoff" in A. 

Figure H] includes numerical evaluations of Ax a for two different sets of assumptions 
including POND. We again see that the functional dependence Axa{E) is roughly the same 
for the different assumption sets. For A = {POND, fixed x 1 }, analytic approximations at 
large E give 

Ax PONDifixed x i = 0.8798E 1/3 - 0.8252 + 0.058E" 1 / 3 + 0.582iT 2/3 - I.eOlir 1 

(39) 

+ 2.25£- 4/3 - 1.16£T 5/3 - 2A1E' 2 + 0(E' 7/Z ) 
To derive (1391) . the first step is to write the relevant energy estimate, (!33l) . as 

E = 7 (— - l) . (40) 
Vl/uv / 

Consider now the combination 

A = 71/uv • ( 41 ) 

The relations (HO!) and ([4"1~]) may be inverted to express 7 and yuv i n terms of E and A. The 
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next step is to use (1351) to express Axpond m terms of E and A. This can be done in closed 
form, but the explicit expression is not very illuminating. The maximum of Ax space ti me is 
attained at a value A = A* determined by the equation 

(-^Ax spacetime (E, A) J =0. (42) 

The asymptotic behavior of A* can be calculated by plugging a large E expansion of the 
form 

A, = Af + \^E-^ + X^E^ 3 + ■■■ (43) 

into (I42pand solving for the coefficients A* term by term. For example, setting the coeffi- 
cients of the first two terms to zero leads to the equations 

-3 V / A^+(l + Al°))cot- 1 V / A^ = 
2(Al°)) 4/3 (l + Al )) + (-2 + Al°)) A« = 0. 

These equations can be solved numerically to give X^ = 0.212 and A* = 0.171. Finally, 
plugging the series expansion (143!) back into our expression for Axpond(-E', A) and expanding 
at large E, one recovers (139]) . 

The method we just described can also be used to obtain the large E behavior of 
Ax spacetime given in ([35]). 



7 Comparing to BDMPS estimates of energy loss 

A standard way of estimating energy loss by hard partons is the BDMPS jet-quenching 
formalism, which is summarized for example in [5j [11] . The energy loss of a hard parton in 
a representation R of the color group SU (N) is 

AE = ±a s C R q(Ax) 2 , (45) 

where Ax is the distance traveled^ Cr is the Casimir denoted C*2(-R) on p. 500ff of [25], so 
that Cf = {N 2 — 1)/2N for a fundamental quark and Ca = N for a gluon. A perturbative 

4 Somc authors prefer to measure distance traveled using a lightconc coordinate, which introduces an 
additional factor of 1/2 into fl5]) . 
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estimate gives 

g pert = ^a 2 iV 2 T 3 . (46) 

7T 

A rule-of-thumb estimate for the strong coupling is a s = 1/2. 

A calculation in M = 4 SYM starting from a Wilson loop definition of q yields [1] 



tew = ^g|^V&^. (47) 

It seems clear that QCD should exhibit a lower value of q because it has fewer degrees of 
freedom: about a third as many as measured by the entropy density. One way to incorpo- 
rating this factor was proposed in [5]: it is to include a proportionality to the square root of 
the entropy in q. Thus 



Scaled LRW ~ y J^j <?LRW ~ 0.63 <? L RW • ( 48 ) 

Because g sca ied still depends on the 't Hooft coupling g YM N, one must fix the value of this 
coupling. An obvious way to do so is to insist that tree-level gluon scattering processes 
should have the same amplitude in Af = 4 gauge theory as in QCD: that is, the tree- level 
couplings coincide, resulting in gy M iV = ® n when a s = 1/2 and N = 3. In quoting a 
numerical value for Scaled lrw in table [□, we have used the "obvious scheme:" g Y M^ = ^ n 
and Tjv=4 = Tqcd- 

An alternative scheme for comparing M = 4 SYM to QCD was proposed in [12]. Instead 
of comparing at fixed temperature, one compares at fixed energy density. This is supposed 
to correct, approximately, for the larger number of degrees of freedom in M = 4 theory, and 
it approximately amounts to setting 3V=4 = T QCD /3 1/4 . Thus, in place of (gH), one would 
have 

^alternative LRW ~ 7^/J<7LRW • (49) 

Also, instead of comparing at fixed tree-level coupling, one chooses in the alternative scheme a 
value (yfy M iV pa 5.5 in order to approximately match the force between a heavy quark and anti- 
quark separated by a distance on order 0.25 fm in a medium at a temperature characteristic 
of RHIC collisions. (Even smaller values of gyu^ can be motivated by comparing SYM and 
QCD at fixed Debye mass [151 •) ^ s shown in tabled], ^alternative lrw ? with <7y M iV = 5.5, is 
essentially indistinguishable from the perturbative estimate fj46|) . 

To compare the BDMPS result to the falling string calculations, we use the operational 
definition (|3J), which amounts to setting the stopping length of a hard gluon equal to the 
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value of Ax that one obtains from (145]) upon setting AE equal to the initial energy of the 
gluon (and, of course, Cr = 3). We further set a s = 1/2 in Q, obtaining 



8E 



(50) 



trough 



3(Ax) 2 ' 



For x in (I50p . we plug in a value estimated from string theory in one of the ways we have 
explained above. Using also ([T|), we find that g roU gti becomes 



We emphatically warn the reader that fl5Tj) is only a rough estimate of the "effective g" 
implied by our falling string picture, because the underlying physical picture is significantly 
different from the BDMPS formalism. Let us review the differences before proceeding to 
extract numbers from (}5~il : 

1. The "gluon" as described by the falling string is off-shell: it follows a time-like trajec- 
tory. This contrasts with the eikonal approximation of lightlike trajectories employed 
in the BDMPS treatment. 

2. In the zero-temperature calculation of [T7], the string worldsheet can be understood 
to arise, in the usual sense of 't Hooft, from a sum over an infinite set of planar 
diagrams contributing to a certain exclusive process; similarly, in our treatment, the 
strong interactions of the gluon with the medium are encoded in the classical dynamics 
of the worldsheet. This again contrasts with BDMPS, which is a partially perturbative 
treatment of radiative energy loss to the medium. It is not obvious to us how to 
translate some aspect of the falling string calculation to the spectrum of radiated 
energy. 

3. We have not included fluctuations in our treatment, so we cannot (yet) give an account 
of the diffusion of transverse momentum similar to the one that is a prominent part of 
the BDMPS formalism. 

Either the obvious scheme (Tqcd — ^V=4 an d 9ym^ = or the alternative scheme 
(Tqcd = 7V=4/3 1 ' 4 and gy M N = 5.5) can be applied to (|5T|) . and we will denote the resulting 
expressions Obvious fail and ^alternative fail- In both schemes we take T QC d = 280 MeV, which is 
a reasonable estimate for central gold-gold collisions at RHIC's top energy, y/s NN = 200 GeV. 
Because x is not exactly proportional to E 1 / 2 , q in equation (|5"T|) depends on E. However, this 




(51) 



21 



obvious scheme 



q (GeV 2 /fm) 

I 



150 



100 



50 



■ • 



■ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 



" ' * ^ • i t H i t i t t t t 4 U 4 i 



LRW 



■ {spacetime, UV} 

▲ {spacetime, trailing, yiR=0.9) 

• {spacetime, fixed x 1 } 

• {POND, trailing, y m =0.9, cutoff} 

• {POND, fixed x 1 } 



10 



20 



25 



E (GeV) 



Figure 6: Evaluations of q as given in ( 151]) under a variety of assumptions, using the obvious 
scheme. The solid blue curve corresponds to the LRW prediction ( 1471) . 



quantity 


Q'pcrt 


<7lrw 


^scaled LRW 


^alternative LRW 


^alternative fall 


Qobvious fall 


value @ 280 MeV 


0.77 


3.6 


2.3 


0.86 


21 


92 



Table 1: Various ways of estimating q. See the main text for details in each one. All values 
for q are quoted in units of GeV 2 /fm. 



dependence is rather weak, amounting at large enough E to a E 1 ^ behavior. Figures [6] and [7] 
show that for 5 GeV < E < 25 GeV, which we take as a representative range of energies for 
hard gluons in the QGP produced at RHIC, q is roughly constant. By fitting x to square- 
root dependence on E for E in the above-mentioned range, and using (1511) . we find 



-{POND, fixed x 1 } _ a GeV 2 -{spacetime, fixed x 1 } _ fi o 

^obvious fall — 110— — <7obvious fall —00 



-{POND, fixed x 1 } _ 9fi 
^alternative fall — ^° 



fm 

GeV 2 
fm 



-{spacetime, fixed x 1 } -, « 

^alternative fall — 10 



GeV 2 

fm 
GeV 2 

fm 



(52) 



We chose the fixed x 1 computation of energy, (1331) . because it avoids the use of IR cutoff and 
thus incorporates fewer assumptions. 

Representative numerical values for the estimates of q discussed here are presented in 
table [TJ The values of g bvious fail and ^alternative fail that appear in this table are the averages 
of the corresponding quantities calculated from POND and spacetime geodesies, as quoted 
in (I52p . We again stress that these averages should be regarded as rough estimates. 

As remarked around (jlj), a recent comparison of parton quenching models to PHENIX 
data favors the range q ~ 7-28 GeV 2 /fm [13]. It is interesting, but in no way conclusive, 
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Figure 7: Evaluations of q as given in ( }5Tl) under a variety of assumptions, using the alter- 
native scheme. The solid blue curve corresponds to the LRW prediction ( 1471) . the dashed 
purple curve corresponds to the "scaled LRW" prediction (HH|) . and the dotted yellow curve 
corresponds to the "alternative LRW" estimate (1491) . 

that the only estimate of q in table [1] that falls within this range is Alternative fail- Of the 
estimates based on [TJ, it is fair to exclude gLRW because it is intended to be for N = 4 
SYM, without any rescalings that would account for the differences between M = 4 and 
QCD. So — with T = 280 Me V — both perturbative estimates and those based on [1] come 
out below the 3a range (J4]), and g bvious fail comes out above. But the T 3 dependence of q 
makes it difficult to pin down the theoretical predictions with much precision. 

We are pleased to see a certain consistency emerging between the falling string calcu- 
lations of this paper and the computations of heavy quark drag and diffusion in j2[ [3l HJ: 
whereas the energy predicted by string theory is too strong when compared in the obvi- 
ous scheme, it is close — though perhaps still a bit high — when compared in the alternative 
scheme |12j . Other comparison analyses have been proposed in the heavy quark setup which 
give pretty good agreement between string predictions and data: see for example [30], in 
which perturbative results are compared between M = 4 SYM and QCD as well as strong 
coupling results. 

The BDMPS approach to energy loss is hardly the only one in common use. Others 
include the higher twist, GLV, and AMY formalisms [3TL |32| |33| 134"] . Each approach makes 
a different set of assumptions. It would be useful to make cross-comparisons with the falling 
string. We leave this task for future work. 
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8 The null string limit 



They key to obtaining an estimate for the stopping length Ax(E) is that the POND trajec- 



tory, which provides an approximate lower bound on Ax, is close to the massless spacetime 
geodesic, which provides an approximate upper bound on Ax: see figure [2] To see why this 
happens, consider the limit v —* 1 with yuv held fixed. A straightforward calculation shows 
that the POND trajectory coincides with the spacetime geodesic in this limit. The key point 
is that the worldsheet becomes a null surface in this limit, so the lightlike tangent vector 
to the POND trajectory is also normal to the worldsheet. In fact, in this limit, k^, and 
n M all coincide up to overall magnitudes. A heuristic way of thinking about this is that in 
the v —* 1, the string is replaced by an ensemble of massless particles, all following critical 
trajectories of the form (126]) (but with different values of K). Signals can't propagate up or 
down the string in this limit: every "bit" of string is causally isolated from every other bit, 
and follows a massless spacetime geodesic. We can develop this heuristic picture by defining 
a scaled version of the string's inverse tension: 



Then, formally, we can take v — > 1 and a' — > oo in such a way that a' scaled remains fixed. 
We will describe this as the "null string" limit, because the string worldsheet becomes a 
null surface. It is a formal limit because when a 1 3> L 2 , stringy corrections to supergravity 
probably become large, so the Ad^-Schwarzschild background is expected to be significantly 
altered. But it captures the key idea that the string tension doesn't matter in the v — > 1 
limit. 

As an application of the null string limit, we can compute the five-dimensional stress 
tensor of the falling string for v — > 1. We are ignoring stringy corrections, so the action we 
start with is 



where Sm is the matter action. The five- dimensional stress tensor t^ u can be defined through 
the equation 



^scaled ^ 



Wl — V 2 , 



(53) 




(54) 




(55) 



so that it enters into the Einstein equation as 



RT - -RG» 



V 




(56) 
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The result will be that the null string's stress tensor is an integral of the stress tensors for 
continuously many massless particles propagating on critical null geodesies. 

First, let's compute the stress tensor of a massless particle, starting from the action ffT9l) 
with m — 0. The result is immediate: 



^.^Ji^-x^)) 1 -—— 
i r dx^ 

d v 6 5 (x fl -X^r]))^— p u (57) 



1 dX» dX u 



V^G J dr) 

1 dX^ 
-.5\x m - X m {z))-—p v 



V^G y v " dz 

where p^ = ~G>iv~p" is the momentum conjugate to X^, and in the last line we have 

specialized to the gauge t] = z. The quantities po and p\ coincide with the expressions in 

x i 

(|23|) when r] = z. Let us now specialize to the critical trajectory (J26J) and set K = + 1, so 

that x l H is the position at which the geodesic crosses the horizon. Introducing the one-form 

^=(-1100 ^p) , (58) 

which is tangent to the particle's trajectory, and employing the gauge rj — z, it is straight- 
forward to show that 

particle = WicleO^) W (59) 

where 



/ n. 1 \ Ezz 

^particle (, i% H ) jj>, 



2 ^-s(t-x 1 H -z H + ^ + £(z)^J 5 (x 1 -xjj-ZH + ^j S(x 2 )S(x 3 ) . 

(60) 

The product of delta functions in (!60|) is simply S A (x m — X m (z)), and we set E = —po- 

Next, we can compute the stress tensor of a string with v < 1 starting from the Nambu- 
Goto action ([6]). The result is 

C mg = ^= / d 2 a 5 5 (x» - X») d a X»P™ (61) 

where is given by ([7]). Taking the null string limit, one finds 

1 1 1 



P» - v^P; = TT-j—jrr "\- h *r h n n * ■ (62) 



null 

string 



2vra^ calcd zjjh\ -h h VT^h, 
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Both rows of (|62|) are proportional to as defined in (158|) . so the matrix has rank 1. 
Observing that coincides, up to an overall factor, with n^, i^, and in the limit v — > 1, 
it already seems inevitable that the stress tensor of the null string will reduce to an ensemble 
of massless particles. Indeed, by plugging fl62|) into flSTj) one finds 

rZn = lim C ing = ~^=0 (x x H - x 1 + z H - ^ ) 6(t - x 1 + ^Mx 2 )6(x 3 )d a X^ 

string f ul1 A/— G V ^ / 

string y 

(63) 

where 

{0 for x < 
(64) 
1 for x > . 

The first delta function in fl63l) enforces the defining relation ffTB"]) of the v — > 1 limit of 
the trailing string. The theta function factor arises because the boundary of the null string 
follows the orbit (JUD with K = |* + 1. Using (152"]) and (1T51), 



1 null 



Tstring^;^)^, (65) 



where 



Tstnn g (^; x^) = - ^ g f x^ - x 1 + ^ - ^ 5(t - x 1 + ^))5(x 2 )5(x 3 ) . (66) 

27r «scaled V Z ) 

To see that the stress tensor of the string is identical to the stress tensor of an ensemble of 
massless particles following critical trajectories, we need only note that 

String (a^jffk) _ 1 L 2 T . , ( rf i. r i \ (a 7 \ 
Ox H 27ra scalcd £,2:^ 

Turning (1671) around, r string is an integral over x\j of r partic i , and since b^b v has no explicit 
dependence on x^, the same relation holds between r s ^ ing and Tp^ rticle . 

A consequence of the discussion of the last few paragraphs is that to compute the con- 
tribution of a falling null string to the expectation (T mn ) of the gauge-theory's stress tensor, 
one can start by doing the analogous calculation for a massless particle and then integrate 
with respect to x\j. The falling null string is, as we have discussed, only a formal approx- 
imation to the finite-tension falling strings of real physical interest. But since its shape is 
known analytically, it seems a worthwhile starting point for an investigation of (T mn ). We 
hope to report on calculations along these lines in future work. 
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So far in this section, we have focused on the limit v — > 1 with ?/uv held fixed. But 
in studying the maximum penetration length of high-energy probes, a different limit is 
appropriate: for fixed energy, we maximize Ax by varying v and yyy so as to hold E 
fixed. For large E, the maximum is attained for 72/uv = ^* ~ ^* > w here = 0.154 for 
A = {spacetime, fixed x 1 } and 0.212 for A = {POND, fixed x 1 }. Thus, the limit of interest 
for estimating penetration length of very hard probes is v — > 1 with 72/uv ne ld fixed — but 
held fixed to different values for spacetime as compared to POND. In this E — > 00, the 
maximum Axpond and the maximum Ax spacc timc do not approach one another: instead, 

AxpoND/AXgpacetimc — > 0.86. 

9 Discussion 

The trailing string of [21 S] is essentially an equilibrium configuration, where energy is lost 
at a constant rate into the plasma, but the quark never slows down because it is infinitely 
massive. A key feature of heavy-ion physics is that many or even most hard partons travel 
only a short distance through the medium before substantially stopping. Our discussion of 
falling strings is a first attempt to incorporate into the trailing string picture the dramatically 
non-equilibrium nature of the physics of energy loss for light partons. It may help the reader's 
intuition to note that in the alternative scheme, x = 2 corresponds to x ~ 0.6 fm, and E = 20 
corresponds to E « 10 GeV, when T = 280 MeV in the real-world plasma. So, according to 
figure HI a 10 GeV gluon stops in a distance of about 0.5 fm. 

Finally, we can consider how the falling string picture might generalize to theories with 
fundamental quarks whose mass is finite. In the construction of [35], fundamentally charged 
quarks come from strings stretched between the D3-branes and D7-branes, where the D3- 
branes create the AdS^, (or Ac/S'5-Schwarzschild) geometry, and the D7-branes are usually 
treated in the probe approximation. At finite temperature, the D7-branes either descend 
to a minimum distance from the black hole horizon if the mass of the corresponding quark 
is sufficiently larger than the temperature, or they extend down into the horizon if the 
corresponding quark is light [3S]- (In any case, provided temperature is constant, the D7- 
branes are static, and apparently stable.) At the risk of oversimplifying, let's ignore the 
geometry of the D7-brane embeddings in ten dimensions and replace them by "flavor-branes" 
that fill AdSs-Schwarzschild either down to some maximal depth < Zjj, or that extend 
across the horizon at z = zjj . The finite- mass trailing string as considered in [2] was assumed 
to end on a flavor brane at the maximum possible depth, i.e. z = z*. A curious property 
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of such strings is that they have a maximum speed, v 



a/1 — z^/zjj, which we alluded to 



following (jSJ). This speed gets smaller as the quark mass gets smaller, and it is in some sense 
zero for quarks that are light enough so that the corresponding D7-brane extends into the 
horizon. It now seems clear that this maximum speed should be understood as the speed 
above which a falling string picture of energy loss must be taken into account. More precisely, 
an energetic quark can be represented as a string coming out of the horizon and extending 
up to a lesser depth, z\jy < z *, than the maximum depth of the flavor brane (if there is one). 
Essentially the same analysis we have given for gluons could be replayed for such strings, 
with the main difference being that the string is no longer doubled, but ends on the flavor 
brane. Admittedly, a fully correct treatment of the D3-D7 construction would involve non- 
trivial motion of the string in the full ten-dimensional geometry, such that the projection to 
AdSs-Schwarzschild would not have the property that the string endpoint travels on a null 
trajectory. However, we are inclined to think that such a motion, which involves the 5*0(6) 
R-symmetry of M = 4 super- Yang-Mills theory, doesn't translate very precisely to QCD. It 
may be that ignoring the ten-dimensional geometry altogether and employing flavor-branes 
in place of bona fide string theory constructions captures approximately the right physics. 
For light quarks, where there is no maximum depth, the only change would then be to replace 
E by E/2: that is, energy loss is half as fast for fundamentally charged light quarks as for 
gluons. This is the same scaling as found in a BDMPS treatment, where according to fj45l) 
the energy scales linearly with Cr, and 



where N = 3 is the number of colors. For heavy quarks, there could be a two-stage process 
of energy loss, where one first has AE cx (Ax) 3 due to falling string dynamics, and then, 
after the string has fallen as far as it can while still remaining attached to the flavor brane, 
one has dE/dx oc p, characteristic of trailing string dynamics. 

We have assumed throughout that string splitting and string joining interactions are 
negligible. In fact they are suppressed by a power of N. In the context of doubled strings, 
the string could split anywhere along its length, and the amplitude to do so is proportional 
to <7 str oc gyMi which is indeed an 0(1/N) effect using 't Hooft scaling. However, this 
suppression may not be enough to make the effect unimportant. The result of one such 
splitting is illustrated in figure Physically, the splitting describes a decoupling of most 
of the hard parton's momentum from its color — a sort of dual hadronization in the fifth 
dimension, where the closed string that carries most of the momentum is a color-singlet 



C F _ N 2 - 1 _ 1 
~C~ A ~ 2N 2 ~ 2 



(68) 
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Figure 8: A doubled falling string which experiences a splitting interaction at a time between 
t — and t = 1. This interaction is suppressed by one power of N. 

glueball in the process of thermalizing. If splitting is significant, it might help the string get 
a little further in the x l direction, but no further than the spacetime geodesic considered 
in section Splitting becomes less easy if the string is not perfectly doubled, as one must 
expect in a more realistic treatment. 

We should keep in mind that falling strings as we have considered them in this paper 
may be considerably idealized in comparison with strings that form as a result of energetic 
collisions. The latter are probably less orderly, and they may be less optimized to travel long 
distances before disappearing into the horizon. This might push our estimates of q upward. 
On the other hand, we have not accounted for fluctuations, nor have we calculated (T mn ). 
Either of these elaborations might be phenomenologically significant, and the latter might 
be affected by the string splitting interactions discussed above. 
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